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Abstract 

Axial vector {J^^ = l"*"^) charmonium and bottomonium hybrid masses are determined via QCD Laplace 
sum-rules. Previous sum-rule studies in this channel did not incorporate the dimension-six gluon condensate, 
which has been shown to be important for 1 and ^ heavy quark hybrids. An updated analysis of axial vector 
charmonium and bottomonium hybrids is presented, including the effects of the dimension-six gluon condensate. 
The axial vector charmonium and bottomonium hybrid masses are predicted to be 5.13 GeV and 11.32 GeV, 
respectively. We discuss the implications of this result for the charmonium- like "XYZ" states and the charmonium 
hybrid multiplet structure observed in recent lattice calculations. 

1 Introduction 

A long-standing problem in hadron spectroscopy is to determine what role is played, if any, by explicit gluonic 
constituents in the hadronic spectrum. Quantum chromodynamics (QCD) suggests the possibility of the existence 
of glueballs which are composed entirely of gluons, as well as hybrids which are composed of a quark, an anti-quark 
and a gluon. An interesting feature of hybrids is that they can have J^^ quantum numbers that are not possible for 
conventional quark mesons. Consequently, the observation of a state with so-called exotic quantum numbers would be 
a "smoking gun" for the existence of hadrons with explicit gluonic content. Hybrids with non-exotic meson quantum 
numbers are possible as well; these could signal their presence through supernumerary states in conventional J^'^ 
channels. In this work we consider the latter scenario. 

Hybrids with non-exotic J^'^ that contain heavy quarks could coexist with conventional heavy quarkonia states. 
The large number of anomalous heavy quarkonium-like states discovered above open flavour thresholds has provided 
an ideal place to look for heavy quark hybrids [3 [H [3J HI [S] . A recent review [5] lists nineteen states discovered 
since 2003, and many of these states are difficult to accommodate as conventional charmonia [7] leading to numerous 
suggestions that some of them may be of an exotic nature (see e.g. Refs. [SlIHlin] for reviews). 

In this paper we use QCD Laplace sum-rules to investigate axial vector (J^'-^ — 1++) charmonium and bot- 
tomonium hybrids. The constituent gluon model [10] was used for the earliest studies of heavy quark hybrids. 
Charmonium hybrids have also been studied using the flux tube model |11| which predicts the lightest charmonium 
hybrids at 4.1-4.2 GeV, as well as lattice QCD [T21[T31[I3] which gives quenched predictions of about 4.0 GeV and 
unquenched predictions of about 4.4 GeV for 1"*"+ hybrid charmonium. Ref. |14| performs a comprehensive study of 
the charmonium spectrum up to approximately 4.5 GeV and finds evidence for a ground state multiplet of hybrids 
which contains the ^ and 1 states, as well as an excited multiplet containing the 1++. As far as we are aware, 
Refs. [ini dH dZ] comprise the only QCD sum-rules studies of axial vector charmonium and bottomonium hybrids. 
Several other channels were examined in this work, and many of the resulting sum-rules exhibited instabilities, lead- 
ing to unreliable mass predictions. The 1++ channel led to well-behaved sum-rules resulting in mass predictions in 
the range 4.7-5.7 GeV for hybrid charmonium and 10.9-11.5 GeV for hybrid bottomonium. 

In recent sum-rule studies of vector (1 ) [18] and pseudoscalar (0 ^) [19] heavy quark hybrids it was shown that 
including the dimension-six gluon condensate can have significant effects on the resulting sum-rules. Specifically, in 
these channels it was found that inclusion of the dimension-six gluon condensate is sufficient to remove the instabilities 
observed in Refs. [H] [161 [17] . With this paper, we explore the effects of the dimension-six gluon condensate on the 
sum-rules for axial vector heavy quark hybrids and provide updated mass predictions. 



In Section [51 we calculate the appropriate two-point function, including leading-order perturbative contributions 
and contributions from the dimension-four and dimension-six gluon condensates. In Section |3l we analyze the sum- 
rules using the single narrow resonance model and then determine ground state mass predictions. Finally, in Sectional 
we discuss the implications of our results for the charmonium-like and bottomonium-like states. With our result for 
the l"*""*" and previous results for the 1 |18| and [19) charmonium hybrid mass predictions, we comment on 
the hybrid multiplet structure identified in Ref. [2]. 



2 Laplace Sum-Rules for Axial Vector Heavy Quark Hybrids 

The axial vector (J^*^ = 1^^) heavy quark hybrids may be examined using the following correlation function [15J 
n^,(g) = J / d^xe^^-^OlT [j,.{x)MO)] |0) (1) 



with Q representing a heavy quark field. Here we examine the transverse part n„ of ([T]), which couples to 1^+ states 

n,.(g) = - 3..) n„(g2) + ^n,(g2) . (3) 

In Ref. [T5l[T^ the perturbative and gluon condensate (aG^) — (aGJ^^G"^) contributions to the imaginary part of 
Hy were calculated to leading order. Here we extend these results by including the contributions of the dimension-six 
gluon condensate (g^G^) = {g^ fabcG^j^^G^^^G'^^^) , which were shown to have important consequences for heavy quark 
hybrid sum- rule studies in different channels [HI [E] . 

First, we verify the leading-order perturbative and (aG^) results [HI [16] for H„. We have opted to calculate the 
full expression for Hi, rather than only its imaginary part as in Ref. |151 116] . This was found to be necessary in order 
to correctly formulate the sum-rules for the pseudoscalar heavy quark hybrid [19 . In addition, verifying existing 
results using a different approach provides a further consistency check of our results. 

The leading-order perturbative contribution to liy is represented in Fig. [TJ We have made use of the Tarcer 
[2Tj implementation of loop-integral recurrence relations and tensor structures [22l [23] to express H„ in terms of a 
small number of basic integrals. Results for these basic integrals are provided in Refs. [24] [SSI [26]. In Z? = 4 + 2e 
dimensions in the MS scheme, the perturbative result is 



SIOOtt^ 
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where terms corresponding to dispersion relation subtraction constants have been omitted. The coupling a and 
quark mass m implicitly depend on the renormalization scale fi in the MS scheme. The generalized hypergeometric 
functions |27j in Q are particularly convenient for sum-rule applications since they clearly reveal the analytic 
structure of H„, namely a branch cut starting at the threshold = Am^. In addition, the imaginary part may be 
easily extracted via analytic continuation of the hypergeometric functions. Doing so, we find 



ImHP"*(g2) = sJY^^l (15 - 35z - 22z^ - 216z^ + 482"*) 

1807r^z"' \ 

. (5) 

+ 15(l-3z + 16z^)log[V^^+^/I] , z>l. 

We find complete agreement between the integral representation for ImHP""^* given in [151 US and (O . 

The leading-order (aG^) contribution to H^, is represented in Fig. [21 Due to the presence of the field strength in 
the current ^ this contribution is most easily calculated using using fixed-point gauge methods (see e.g. Ref. [2H] 
for examples of this technique). However, plane wave methods could be also be used as they have been proven to be 
equivalent to fixed-point gauge when gauge-invariant currents such as ^ are used [29j . For the {aC^) contribution 
we find 

n?^(g') = + 2^)2^1 (1' 1; 5/2; z) , (6) 
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Figure 1: Feynman diagram for the leading-order perturbative contribution to n„. The current 
is represented by the (8 symbol. This and all subsequent Feynman diagrams were created with 
JaxoDraw [20] , 



where non-physical terms corresponding to dispersion relation subtraction constants have been omitted. The imagi- 
nary part of ([6]) is 



tGG/ 2^ m2(aG2) , y/z-l 



Imn:r^(g^)- (1 + 2z) , z>l (7) 

This again agrees with the result given in Refs. [151 116) . 



Figure 2: Feynman diagram for the leading-order (qG^) contribution to n„. 

Finally we consider the dimension-six gluon condensate contributions which were not calculated in Refs. [151 116) . 
These are represented by the diagrams in Fig. [31 Again utilizing fixed-point gauge methods, we find 



" 1^2 [2^(2 - 9^ + 6^') - 4z (z - 1) (3z ~ 1)] 2^1 (1, 1; 5/2; z) 

iiO/TT {Z — ij 

+ 1152^2(,li)2 [3(17^ - 9)(z - 1) - 3(17 - 46z -I- 21 z^)] . 
The resulting imaginary part of ([8]) is 



Imn«°«(g2) = 38i4^!\)2 [2(1 - 3z)(z - 1) + (2 - 9z 6z^)] , z > 1 . (9) 

The singularity at z = 1 in © must be dealt with carefully. Although ^ can be extracted from the purely 
hypergeometric terms in it is not well-defined at z = 1. This problem is addressed through inclusion of the non- 
hypergeometric terms in (jSl). These contribute compensating terms which ensure that the contributions of Il'^'^'^{q^) 
to the sum-rules are well-defined when Imn^'^^(g^) is integrated from z = 1. Thus the imaginary part ([9]) by itself 
is insufficient to construct the contribution of the dimension-six gluon condensate to the QCD Laplace sum-rules. 

Now that we have calculated the correlation function, we can proceed with the QCD Laplace sum-rules analy- 
sis fSnilSI] (for reviews of the methodology see e.g. Refs. [311 [33])- Utilizing the standard resonance plus continuum 
model for the hadronic spectral function, the Laplace sum-rules take the form 

1 1"^ 

£^^° (r, so) = - / t"" exp [-tr] p^^^it) dt , (10) 

where is the hadronic threshold. The quantity on the left hand side of (flUl) is given by 

1 r n 1 t°° 

£^^° (r, so) = -B {^if Q'^^n, (Q2) - - / t'' exp [-tr] lmR,{t) dt (11) 
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Figure 3: Feynman diagram for the leading-order {g^G^) contribution to Hi,. Additional diagrams 
related by symmetry are not shown. 

where sq is the continuum threshold, = —q^ is the Euclidean momentum and 11^ the axial vector 

heavy quark hybrid correlation function. The Borel transform operator B is closely related to the inverse Laplace 
transform [34] 

^-B[f{Q')]^F{r)=C-'[f{Q')] , (12) 

b-j-ioo 
b—icc 

where b is chosen such that / (Q^) is analytic to the right of the integration contour in the complex plane0 Therefore 
any terms in the full expression for the correlation function n^, (Q^) that contribute to the inverse Laplace transform 
(|13p must be included in the construction of the Laplace sum-rules. The singular terms in ([8]) that do not contribute 
to the imaginary part (j9|) fall into this category, and thus they are an essential element of the QCD Laplace sum-rules. 
Using the results for the leading order perturbative ([5]), (aG^) Q and (g^G^) ([8]), ([9]) contributions, we find 

£^^° (r, so) = / [Imnf'-* {'im'^x) + ImUf^ {4m^x)] exp (-4mVx) dx 

Ji , . 

/ Imn^^^(4TO'^x)exp (-4m^Ta;) dxH ^exp(-4m^r) 



lim 



Cf^^(^r,so) = ~^Cr''{r,so). (15) 

The terms involving rj in (|15p render the integration in ()14p well-defined for the x 1 (77 — > 0) limit, and are natural 
consequence of the inverse Laplace transform of the full expression Again, we stress that this expression cannot be 
obtained with Imlljf '-^'-^ alone. As before, the mass and coupling in ([T4|) and ([T5|) are functions of the renormalization 
scale /i in the MS-scheme. After evaluating the t derivative in (|15p . renormalization group improvement may be 
implemented by setting /i = l/^/r [36] . 

3 Analysis: Mass Predictions for Axial Vector Heavy Quark Hybrids 

In order to extract ground state mass predictions for the l'^^ heavy quark hybrids, we use a single narrow resonance 
model 

^p^^'^it) = fS {t - m2) . (16) 



Using this in Eqn. (jlOp yields 

>C?''°(r,5o) = /'M2'=exp(-MV), (17) 
from which the ground state mass M can be determined via the ratio 

= . (18) 



'^Ref. |35| contains detailed examples applying inverse Laplace transform techniques in sum-rule calculations. 
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It should be noted that the narrow resonance model (fTC|) results in a larger mass prediction M compared to resonance 
models including width effects [37]. Additionally, an upper bound on the ground state mass prediction M can be 
obtained by taking the limit as sq — >■ oo. The resulting upper bound on M is quite robust as it does not depend on 
the resonance model or how the QCD continuum is modelled. 

To extract a ground state mass prediction the sum-rule parameters must be fixed. In the interest of self- 
consistency, we have chosen to utilize sum-rule estimates of quark masses. For the charm and bottom quark masses 
we take 

TOc(/i = = mc = (1.28 ±0.02) GeV, (19) 
TOf,(^ = mb) =mh = (4.17 ±0.02) GeV , (20) 

corresponding to the full range of MS charm and bottom quark mass estimates of Refs. [38l [39l HQI [41] and in 
agreement with the ranges recommended by the Particle Data Group 42 . We have used one- loop MS expressions 
for the coupling and quark masses. The coupling is evolved from the t and Z mass for charmonium and bottomonium 
hybrids, respectively: 

= — TaIT^^^ / . X ' " = 0-33 ; (21) 



127r 6 \ M 



The numerical values of a [Mr) and a {Mz) are taken from |43|, and we use Particle Data Group values of the r and 
Z masses [IJj. At one- loop order, the MS charm and bottom quark masses are given by 



_ / a{n) 



12/25 

TOc(Ai) = rUc [ ) , (23) 

/ ( \ \ 12/23 

m,{^,) = mA^^] . (24) 

For the purposes of the sum-rule analysis we set fi — 1/ ^ as described above. We use the following values of the 
QCD condensates, extracted from heavy-quark systems [40] : 

(.g^G^) = (8.2 ± 1.0) GeV2(aG'2) (25) 
(aG^) = (7.5 ± 2.0) x 10"^ GeV^ . (26) 

We find that 77 = lO^"' is sufficient to numerically evaluate the limit in ([T4|) . 

Now that the numerical values of the physical parameters have been fixed, we may proceed with the sum-rule 
analysis beginning with hybrid charmonium. First we must establish a region of validity for the sum-rule analysis. 
To do so, we follow Ref. 131 and define the functions 



£Q^"(r,.o)/£Q^°(r,.o) 

/cent (r, .0) - ^QCD ^) /^QCD (27) 



/:Q^"(r,.o)/£?^"(r,.o) 

;pow(r,.o)- cFr{r,s,)llir\TM ' ^ 

where >C^"* represents the perturbative contributions to ([T?)) and ([15]). The functions ([?7)) and (|^5]) measure the 
relative importance of the respective continuum and non-perturbative contributions to the sum-rule. These two 
functions can be used to constrain the Borel parameter r and define a window of reliability for the sum-rule. 
Inspired by Ref. [31], we impose the constraints /cont > O-^ (i.e., continuum contributions must be less than 30%) 
and |/pow~ 1| < 0.15 (i.e., non-perturbative contributions do not exceed 15%). The purpose of these constraints is to 
control uncertainties associated with non-perturbative effects (such as truncation of the operator-product expansion 
at dimension-six) and the continuum model. Figure [J] depicts the functions /cont and /pow for the optimal value of 
So which is determined below. With these criteria we find the region of validity to be 5.3 GeV^ < l/r < 7.3 GeV^. 
We have also performed the analyis in the pole scheme, with a charm quark pole mass of vnP°^^ = 1.71 GeV 42 . As 
in the pseudoscalar charmonium hybrid analysis [H], we find that the sum- rule window closes rapidly in the pole 
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Figure 4: The quantities /cont (t, so) (solid line) and fpov (t, sq) (dashed line) for hybrid char- 
monium are shown as a function of the Borel scale 1/r for the optimized value sq = 33.0 GeV^. 
Central values of the QCD parameters have been used. 

scheme, and hence the MS scheme is more suitable for this analysis. The advantage of the MS scheme is also seen 
in Ref. [H]. 

Once the sum-rule window has been identified, the optimal sq value is determined as follows. First, the lowest 
value of So where the mass prediction (|18p stabilizes (exhibits a minimum) within the sum- rule window 5.3 GeV^ < 
l/r < 7.3GeV^ is identified. In this case, this criteria establishes the constraint sq > 32GeV^. Figure [5] shows the 
mass prediction ()18|) within the sum-rule window for several values sq. Second, we define 

and then search for the value of sq that minimizes p9|) in the region where sq > 32 GeV^ and 5.3 GeV^ < l/r < 
7.3 GeV^. This results in an optimal sq = 33GeV^ and a corresponding charmonium hybrid mass prediction of 
5.13 GeV. Note that the limit as sq — > oo cannot be used here to obtain a demonstrable upper bound on the 
charmonium hybrid mass since the mass prediction (1181) does not stabilize within the sum-rule window, as can be 
seen from Fig. [5l 

-^Xl/Xo (GeV) 
6.2 - 
6.0 - 

5.8 - ^ - - ' " 

5.6 ~ 
5.4 - 
5.2 - 
5.0 - _ 

48 t • \ • • • • \ • • • • \ • • • • \ • • . 1/t (GeV-) 

5.5 6.0 6.5 7.0 

Figure 5: The ratio C^'^^ (t, sq) /£^'~^'^ (t", sq) for hybrid charmonium is shown as a function of 
the Borel scale l/r for the optimized value sq = 33 GeV^ (solid curve). The ratio is also shown for 
So = 38GeV^ (upper dotted curve), so — 28GeV^ (lower dotted curve) and so — > oo (uppermost 
dashed curve). Central values of the QCD parameters have been used. 

We now estimate the uncertainty in the charmonium hybrid mass prediction due to uncertainties in the QCD 
input parameters. Interestingly, the uncertainty in the mass prediction is dominated by (aC^) (|26p . while the 
uncertainties due to the charm quark mass (|19p and {g'^G^) (|25p are significantly smaller. This is in contrast with 
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the pseudoscalar heavy quark hybrid, where the uncertainty in the mass prediction is dominated by the dimension- 
six gluon condensate. We have made no attempt to estimate contributions to these uncertainties from higher loop 
effects. Adding the uncertainties in quadrature, we predict the charmonium hybrid mass to be 5.13 ± 0.25 GeV. 
This prediction is in good agreement with the range of results of 4.7-5.7 GeV found in Refs.[Tni [T^, all of which 
were derived from stable sum-rules that did not include effects of the dimension-six gluon condensate. Thus we can 
conclude that the dimension-six condensate is not as significant in the 1"*"+ channel as it is in the ^ and 1 
channels of hybrid charmonium. 

The sum-rule analysis of hybrid bottomonium is very similar. The sum-rule is reliable in the region sq > 145 GeV^ 
and 7.8GeV^ < l/r < 25.0GeV^. The functions /cont (l?7l) and /pow are shown in this region in FigureHl and 
we again use the constraints 0.85 < /pow < 1-15 and /cont > 0.7. As in the hybrid charmonium analysis, it is not 
possible to obtain a demonstrable upper bound on the mass prediction since the ratio (jl8p for sq — > oo does not 
stabilize within the sum rule window. Using (|29|) . we find the optimal sq — 150 GeV^. Figure [7] shows the mass 
prediction (|T8| within the sum-rule window for several values of sq. The uncertainty analysis again shows that the 
error in {aG'^) dominates the error in the bottomonium hybrid mass prediction. Adding the errors in quadrature, 
we predict a bottomonium hybrid mass of 11.32 ± 0.32 GeV. This result is in good agreement with the range of 
10.9-11.5 GeV for hybrid bottomonium predicted in Refs. [T51 117) . 




Figure 6: The quantities /cont (t, so) (solid line) and /pow (t, so) (dashed line) for hybrid bot- 
tomonium are shown as a function of the Borel scale l/r for the optimized value sq — 150 GeV^. 
Central values of the QCD parameters have been used. 
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Figure 7: The ratio Ci (t, so) /^o (t, so) for hybrid bottomonium is shown as a function of 
the Borel scale l/r for the optimized value so — 150 GeV^ (solid curve). For comparison the ratio 
is also shown for so = 170 GeV'^ (upper dotted curve), so = 130 GeV^ (lower dotted curve) and 
So — ^ oo (uppermost dashed curve). Central values of the QCD parameters have been used. 
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4 Conclusions 



In this paper we have studied axial vector (J^*^ = 1++) heavy quark hybrids via QCD Laplace sum-rules. We 
have calculated the full expressions for the leading order perturbative and (aG^) contributions to the correlation 
function, and have noted that the corresponding imaginary parts of these expressions agree with the results given in 
Refs. [121 [TT]. For the first time, we have also determined the contributions from (g^G^), which were not included in 
previous work. For these it was shown that the imaginary part alone was insufficient to formulate the QCD Laplace 
sum-rules, and the full expression for the (g^G^) contribution to the correlation function was needed. 

In Refs. [151 [17], many of the sum-rules for various J^'~'' heavy quark hybrids exhibited instabilities, and hence 
the resulting mass predictions for those channels are unreliable. Recent sum-rule analyses of vector (1 ) [TB] and 
pseudoscalar (0 '") (19j heavy quark hybrids have shown that the inclusion of the dimension-six gluon condensate 
stabilizes the sum-rules for these channels. Although no instabilities were found for 1++ heavy quark hybrids in 
Refs. [ISl [U], it is nevertheless of interest to examine the effects of the dimension-six gluon condensate in this 
channel given the significant effect observed in other channels. Including the {g^G'^) contributions in our analysis 
results in the predictions of 5.13 ± 0.25 GeV for hybrid charmonium and 11.32 ± 0.32 GeV for hybrid bottomonium. 
Our results are in agreement with the predictions of Refs. [HI [13 which ranged from 4.7-5.7 GeV and 10.9-11.5 GeV 
for axial vector charmonium and bottomonium hybrids, respectively. The uncertainties in the mass predictions are 
dominated by the uncertainty in (aG^) while the uncertainty due to (g^G^) is less significant, in contrast to the 
uncertainty in the pseudoscalar hybrid mass predictions [19) . This, together with our agreement with the mass 
predictions of Refs. [HI [17], suggests that unlike the vector [18] and pseudoscalar [19^ channels, the effects of the 
dimension-six gluon condensate are less significant for the axial vector channel. 

To date, all of the charmonium-like "XYZ" states have been discovered in the mass range 3.8-4.7 GeV jB]. Clearly 
our prediction of 5.13 GeV does not support the identification of any of these states as an axial vector charmonium 
hybrid. The first discovered charmonium-like state was the X(3872), whose possible J^'~' assignments are 1++ or 
2 ^ [45l[46], although the 1++ option is strongly favoured [9]. Many different proposals have been made regarding 
the nature of the X(3872): a conventional charmonium state, a D^* molecule, a tetraquark and a hybrid (see 
e.g. Ref. [S] for a review). The hybrid interpretation was suggested in Ref. [JT], but has now been largely set aside. 
The reason for this seems to be that both flux-tube model 11 and lattice QCD [T21[I31[T1] predict that the lightest 
charmonium hybrids have masses significantly greater than that of the X(3872). If its quantum numbers are shown 
to be 1++, our mass prediction of 5.13 GeV is in agreement with the results of other theoretical approaches that 
seem to preclude a charmonium hybrid interpretation of the X(3872). 

It is interesting to note the large difference between the predicted masses of 3.82 GeV for ^ and 4.12-4.79 
for 1 [18] hybrid charmonium compared to the 1"*"+ prediction of 5.13 GeV. In Ref. [14] it is suggested that ^ 
and 1 are members of a ground state charmonium hybrid multiplet, while 1++ is a member of a multiplet of 
excited charmonium hybrids. Although the mass splittings are significantly larger, the present result and those of 
Refs. |181 119] seem to be in approximate agreement with this multiplet structure. Future work to update remaining 
unstable sum-rule channels in Refs. [151 116[ [IT to include the effects of the dimension-six gluon condensate would 
clarify the predictions for the spectrum of charmonium hybrids from a QCD sum-rules standpoint. 
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